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Abstract. We propose a method for quantum state transfer in spin chains using an 
adiabatic passage technique. Modifying even and odd nearest-neighbour couplings in 
time allows to achieve transfer fidelities arbitrarily close to one, without the need for 
a precise control of coupling strengths and timing. We study in detail transfer by 
adiabatic passage in a spin-1 chain governed by a generalized Heisenberg Hamiltonian. 
We consider optimization of the transfer process applying optimal control techniques. 
We discuss a realistic experimental implementation using cold atomic gases confined 
in deep optical lattices. 
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1. Introduction 

Reliable quantum communication between distant nodes in a quantum network is of 
uttermost importance for quantum computation and information [1]. For relatively 
short distances, since the pioneering works of S. Bose [2] and V. Subrahmanyam [3] 
spin chains are considered to be good candidates to perform this task. However, due 
to dispersion in the free evolution, the fidelity of the transmission through such a chain 
decreases as the number of spins is increased. This problem was first circumvented 
through local engineering of the (nearest-neighbor) couplings in the chain to obtain a 
perfect channel [1]. Since then, many schemes have been proposed to maximize transfer 
fidelity and practicability. Coupling the outer spins only weakly to the chain allows 
to perturbatively obtain an effective Hamiltonian for the two outer sites only [5l [6]. 
Also, dispersion can be minimized by suitably encoding a quantum state in several 
neighbouring spins [3 E]- If the quantum state is encoded in two parallel, uncoupled 
spin chains ('dual-rail' encoding), then a suitable measurement at the end of the chain 
allows to conclude that the state has been transferred successfully [9]. All these schemes 
rely on chains or ladders of spin-l/2s. Transfer of gaussian states has been studied in 
chains of harmonic oscillators [IDl [H] • For most protocols, no temporal control of the 
couplings inside the chain is necessary. However, usually it is inevitable to control the 
first and the last coupling in order to start the transport and to receive the state at the 
other side, respectively. Especially the control of the last coupling often is crucial to 
retrieve the state with maximal fidelity. 

Here we study quantum state transfer in one-dimensional (ID) spin-1 systems, requiring 
a (limited) temporal control of all the couplings in the chain. We focus on the most 
general isotropic spin-1 Hamiltonian, the generalized Heisenberg model. Its rich phase 
diagram, which includes several paradigmatic anti-ferromagnetic phases of strongly 
correlated many-body systems, has been intensively studied for the last two decades 
(see, e.g., [12] and references therein). Still, properties of ground states and low 
excitations are under debate in certain regions [I2]. Recently, this type of Hamiltonians 
has attracted a renewed interest because their clean realization is possible by loading 
spin-1 atoms (e.g., ®^Rb or ^^Na) into a deep ID optical lattice [131 [HI [15]. Quantum 
state transfer through spin-1 chains governed by the generalized Heisenberg Hamiltonian 
has been studied in [16] for the case when the chain with all couplings fixed and 
identical is prepared in its ground state. Like in the spin-1/2 case, the free evolution 
after a single spin is coupled to the chain presents dispersion. Therefore the transfer 
is generically non-perfect. However, certain (anti-ferromagnetic) regions in the phase 
diagram give rise to particularly high transport fidelities. The transfer mechanisms in 
these cases are strikingly different from those in spin-1/2 chains which usually start from 
a ferromagnetic state. 

Here we propose a scheme that uses adiabatic passage techniques to obtain arbitrarily 
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perfect state transfer 0. It requires the ability to modify even and odd couplings in the 
chain independently. Such a control can be realized for atoms confined in an optical 
lattice via super-lattice techniques. By choosing the duration of the process long enough, 
the scheme allows to increase the transport fidelity arbitrarily close to unity, without the 
need for a precise control of the couplings. A transfer scheme using adiabatic passage 
techniques is inevitably slow as compared to the speed of propagation of a low excitation 
in the system with all couplings active and identical (the same happens for schemes using 
local engineering of couplings or working in the perturbative regime). We will, however, 
demonstrate that optimal control techniques allow to significantly reduce the duration 
of the transport process. 

The paper is organized as follows: In Section [2] we introduce the spin-1 generalized 
Heisenberg Hamiltonian with bilinear-biquadratic nearest-neighbor couplings. In 
Section [3] we provide a description of the transfer scheme and discuss in detail its 
performance. In Section |4] we consider optimization via optimal control techniques. A 
realistic implementation of the proposed scheme in a system of ultracold atoms confined 
in a deep optical lattice is discussed in Section [5l Finally, we conclude in Section [6l 



2. Generalized Heisenberg Hamiltonian 



We consider a one-dimensional chain of spin-Is, coupled through the most 
general isotropic ID Hamiltonian (the generalized Heisenberg, or bilinear-biquadratic, 
Hamiltonian) with nearest-neighbor interactions 

Hbh = <yy . cos6 (ji J i+i) + sine (ji J. i+i) (1) 
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Here, Jj denotes the vector of spin-1 operators for site i and a > fixes the coupling 
strength. The first term in f^bb is the usual linear Heisenberg interaction. The second, 
quadratic, Heisenberg term does not appear in spin-1/2 systems, as in that case any 
power of (JjJj+i) can be expressed as a combination of a constant and a linear term. 
The properties of ffbb often become more intuitive if it is expressed in terms of two-site 
operators Pa''^~^^\ projecting sites i, i + 1 onto manifolds with total spin S. Using the 



relation 



we obtain 
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I During completion of this manuscript, an article of T. Ohshima et al. has appeared [29] which 
discusses quantum state transfer in spin-1/2 chains through adiabatic dark state passage. 
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with Xq = —2 cos 6 + A sin ^, Ai = — cos 6 + sin 6, and A2 = cos 6 + sin 6 [13]. Note that 
is always possible to set the smallest Aj to zero by adding a multiple of the identity to 
Hhh- 

Properties of this spin-1 Hamiltonian, Eq. ([T]), especially of its ground-states, have 
been extensively studied (see, e.g, [121 [131 [IS [HI [HI [19]. It is the interplay between 
the bilinear and biquadratic coupling which gives rise to a rich manifold of ground 
states. The ground state is ferromagnetic {i.e., magnetized) for ^ti < 9 < and 
anti-ferromagnetic {i.e., non-magnetized) otherwise. Here we concentrate on the fully 
dimerized phase, 9 = —it/ 2. To better understand its properties, let us start by 
considering a system of only two coupled sites. Since for 9 = —it/2 the two-site 
Hamiltonian takes (after a constant shift) the simple form -^bb*"*"^^ ~ —Pq'^^^\ the 
energy is minimized by a singlet state |s)i2 = (|l)i| — 1)2 + I — l)i|l)2 — |0)i|0)2) /V^. 
Here we denote the three eigenstates of as | ± 1), |0). For more than two sites, such 
a configuration cannot be repeated on neighboring bonds. Still, for an even number of 
sites and open boundary conditions, a good caricature of the ground state is given by 
a dimer (or valence bond) state, which has a singlet on each second bond. Though we 
concentrate here on 6* = —it/2, the scheme we propose works in a region of the phase 
diagram around this point. 

3. Transfer scheme 

To achieve the state transfer, we assume an odd number of sites, = 2n + 1, and 
modify the Hamiltonian Eq. ([ID as follows (here and in the following we set 9 = —tt/2): 



i.e., we introduce the parameter x G [—1,1] to be able to adjust the ratio between 
even and odd couplings. For x = 1 (a; = — 1) only even (odd) couplings are non-zero. 
Realization of such a Hamiltonian for ultracold atoms (as ^^Na) confined in an optical 
lattice is discussed in Sec. [H 

Let us start by analyzing in detail a system of = 3 sites. Setting initially x = 1, only 
the coupling between sites 2 and 3 is turned on, while the first site is decoupled. The 
(degenerate) ground state then consists of an arbitrarily oriented spin at site 1 and the 
two other spins paired into a singlet. Thus 



On the other hand, if x = —1, only the first two spins are coupled. In the ground state 
they then form a singlet, while the last spin is in an arbitrary state 




N-l 



(4) 



1=1 



(5) 



I^L-l(0')) = I^)l2|0')3. 



(6) 
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We will now show that \ip^^i{(f))) and \ip^^_i{4>)) are connected via an adiabatic path 
as X is changed continuously from 1 to —1. Then, a protocol to transfer an arbitrary 
state 10) between the two spins at the end of the chain can be constructed as follows: 
(i) at t = tstart choose x = 1 and prepare the system in its ground state, with the first 
spin initialized to the state |0) to be transported: \ipit = tstart)) = |^S=i(0))i (ii) change 
X to reach x = —1 at t = tfinai- If the total duration of the process T = tgnai — ^start 
is long enough to have a sufficiently slow change of x, then according to the adiabatic 
theorem [20] the system remains in the same band of instantaneous eigenstates, and 
li^if = ^finai)) ~ |'0°=_i(0))- Note that a precise timing for decoupling and/or reading 
out the final spin after the transport is not crucial here (as it is for transport schemes 
using fixed couplings [21 HI [3]), as the coupling of the last spin is adiabatically set to 
zero. 

We take |0) = = |1) and denote \L) = |l)i|s)23, |-R) = |s)i2|l)3. Observing 
Pq^'^^\L) = \R)/S, we find that Hi,\y{x) only acts in the subspace spanned by {\L), \R)}. 
As {L\R) = 1/3, we introduce an orthonormal basis as 

{\b,) = + \R)), \b,) ^ y^{\L) - \R))}. (7) 

The Hamiltonian, restricted to this basis and shifted by —al/2, reads 



Corresponding eigenvectors are 

1-).= cos^|&i) + sin^|62), (9) 

1+).= sin cos ^162), (10) 

with tan,^(x) = 1\plx. The energy eigenvalues are e=p = =paVl + 8x^/6. We have 
= \L) and |— )^^_i = |-R). Since due to the adiabatic theorem the system 
remains in the same band of instantaneous eigenstates when the parameter x is changed 
slowly from +1 to —1, population can be transferred from \L) to \R) through adiabatic 
passage. Notice that the time dependence in li)^; and e± is explicit through the control 
parameter x(t). 

A characterization of the adiabatic regime necessary to obtain transport with high 
fidehty is provided by the condition [21] 
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Figure 1. Excited state population p^^ obtained from Eq. (fT4|) as a function of the 
duration T of the process. 



Here it is assumed that x is changed at a constant rate, i.e., dx/dt = 2T ^, and T then 
is the total duration of the process. Furthermore, 

cW..H-i|_>^ = _^ (12) 

is the couphng between \—)x and \+)x, and 

A(a;)^e,-e_ = ^v^rT8^ (13) 



is the energy difference of the two levels. Using Eqns. ( 1121) and ( 11311 we obtain 
T ^ 3^/2 /a. In this limit, first order perturbation theory allows to obtain an estimation 
of the final population p+ of |+) due to the non-adiabatic coupling [22] {h = 1): 



/ dx exp [ iT dx' A{x') j c{x) 

Jx=l \ Jx'=l / 



(14) 



The excited state population p^ is plotted in Fig. [T] as a function of the total process 
duration T. For T large, the term exp(T J^,^-^^dx' A{x')) oscillates rapidly, and due to 
destructive interference the excited state population nearly averages to zero after one 
cycle. For this reason, p+ generically decreases as T is increased. The additional dips 
in the curve correspond to values of T for which interference is such that p+ is reduced 
after the whole process. 

For 10) = |0), I — 1) the same effective Hamiltonian (Eq. [HI in the respective basis) is 
obtained. Thus any state is transferred with the same efficiency and the same dynamical 
phase as |1). This follows from the rotational symmetry of the Hamiltonian as well as 
of the initial state of the chain. Additional terms in the Hamiltonian breaking the 
rotational symmetry will in general introduce differences in the dynamical phases and 
a change in the adiabaticity condition. As an example, let us discuss the presence of a 
global magnetic field term oth^^ in the Hamiltonian (j4j). We assume h to be small 
in order to still have a dimerized ground state. The energy of the states \^)x then is 
shifted by ahStot, where — 1 < ^tot < 1 is the total magnetization of the chain. Both, 
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\+)x and \—)x, are shifted by the same amount, such that the couphngs and thus the 
adiabaticity condition are not altered. However, the shift depends on the initial state 
10), and a differential dynamical phase is introduced. This can be corrected after the 
transfer by applying a local unitary to the last spin. Dynamical phases from fluctuations 
of h around a constant value nearly average out, as long as they happen on a time scale 
short compared to the total process duration. Large field fluctuations on short time 
scales additionally might spoil the adiabaticity, leading to larger couplings to excited 
states. 

We can immediately generalize this scheme for a larger chain with an odd number of 
spins. For x = 1, only even couplings are turned on. Then the ground state consists of 
a decoupled spin at site ^ = 1, and all other spins are paired into singlets: 

IV^Ll(0)) = \(P)\s)2,3 ■ ■ ■ \s)2i,2i+l ■ ■ ■ |s)2n,2n+l. (15) 

On the other hand, if x = — 1, such that only odd coupling are active, then the ground 
state has again n pairs of singlets, now with the last spin decoupled: 

IV'L-ll^')) = \s)l,2 ■ ■ ■ \s)2i-l,2i ■ ■ ■ |s)2n-l,2n|0')- (16) 

The magnetization M = Jj) of \ipx=-i) and \ip^^_^i) is completely fixed by the first 
and the last spin, respectively, as the rest of the chain is non-magnetized. Also, the 
magnetization is conserved under the rotationally invariant set of Hamiltonians Hiyi^{x), 
and moreover it can be shown that for any — 1 < x < 1 the ground state is unique 
for a given magnetization [18]. In a finite system, it is furthermore separated by a 
gap from the first excited state. We can thus apply the same protocol as for = 3 
for transporting an unknown state of the first spin to the last one. An example of the 
population transfer via adiabatic passage is depicted in Fig. [2] for a chain of 5 sites. The 
first spin is initialized to |0) = |1). Populations of |1) for site i, Pi{t) = i{l\pi\l)i {pi is 
the reduced density matrix of site i), are shown as a function of time as x is changed 




Figure 2. The populations Pi{t) = i{l\pi{t)\l)i for a chain of 5 sites as a function of 
time, changing x from 1 to —1 at constant rate. 
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Figures. Logarithmic plot of the error e = 1 — F= 1 — 7v(l|pw(i)|l)Ar versus velocity 
T/{N — 1) for various lengths of the chain {6 = —tt/2). The couplings are changed 
linearly: x{t) = (T — 2t)/T . Time evolution has been performed using matrix product 
state algorithms j25| . 



linearly from 1 to —1 in a time interval T = 80a~^ = !)• As expected, a significant 
increase of the population of |1) is only present for odd sites. For perfect adiabatic 
evolution, we should find pi{t) = pN~i{T — t). The deviations (wiggles for t > T/2) 
stem from non-adiabatic transitions. 

Let us discuss the conditions required to be in the adiabatic regime (Eq. ( fTTl) ) for a 
chain of length N. Using \{ij^{x)\d/dx\ij^{x))\ = |(?/;^(x)|(rf^bb/c^a;)|V^°(x))/A(x)| and 

together with | I-Pq'*'*'^! | = 1, c(x) can be bound from above as c{x) < {N — 1)/A(x). 
The minimal gap min^; A[x) is more difficult to estimate. For odd and 6 = —ixjl it is 
know from the equivalence of tih\,{p^ = 0) to an XXZ model, that in the thermodynamic 
limit, i.e., for N — > oo, the gap vanishes: lim7v-»oo Amin = [23]. For small chains, we 
find that the min^, A{x) decreases approximately linear in 1/N. Then the minimal value 
T necessary for being in the adiabatic regime grows approximately with the third power 
of the number of sites. 

The error of the transfer process is given by e = 1 — F, with F = i\i{l\pi\f{T)\l)N being 
the transfer fidelity. In Fig. [3l the error e is plotted for different durations of the process 
and for chains of various lengths. Clearly, e decreases as the transfer process is made 
slower. The maximal velocity T/ (A^ — 1) permitted for state transfer with fixed maximal 
error decreases as N is increased. 
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4. Optimization of the transfer fidelity 



For the Hamihonian considered here, the decrease of the transfer velocity as the number 
of sites is increased cannot be avoided, since the gap to the first excited state closes for 
N —>■ oo. For small chains {N < 11) we have performed an exact diagonalization of 
Hhh{x) to obtain c(x) and A(x). Figs. H] (a,b) show typical values of c{x) and A(x) for 
the coupling of the ground state to the first excited state (both with the same total spin). 
The gap A(x) decreases approximately linear for |x| varying from 1 to 0, while c(x) is 
well described by a Lorentzian function. Clearly, both are sharply peaked at a; = 0. 
This suggest to use optimal control techniques to increase the (mean) transfer velocity 
by adapting the change of the spin-spin couplings, i.e., dx(t)/dt, to the instantaneous 
energy difference A{x) and the coupling c{x). In [22], a method has been developed 
to obtain an optimized path x = Xopt{t) which allows to suppress excitations to other 
states, i.e., to reduce the transition probability 



p{T) 




dt'Aixit'))^ cixit))^ 



To this aim, we replace J^,^^dt'A{x'{t)) = ^T{t) (where Tq is chosen to have r(0) = 
and t(T) = 1), such that the integral becomes the Fourier transform of a function 
n(r) = ^^^^c(x(r)). Now, for n(r) we choose a Blackman pulse [23] to reduce the 
side lobes of the Fourier transform. The shape x{t) is then obtained from solving two 
differential equations in order to obtain first x(r) (from inverting the equation for u{t)) 
and subsequently r(t). 

An example of a path x(t) obtained in such a way is plotted in Fig.|4](c). The dependence 
of the transfer error on the mean velocity T/{N — 1) if such an optimized path is used, is 
plotted in Fig. [51 Compared to the results for a linear change of x in time, displayed in 
Fig. [31 the error for a given mean velocity is typically more than one magnitude smaller 
for the optimized path. 




Figure 4. (a) Coupling c{x) and (b) gap A(x) to the first excited state as a function 
of the parameter x for 6 = —tt/2 and N — 7 sites, (c) The optimized path x = Xopt{t), 
obtained from the curves of (a,b). 
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Figure 5. Error e = 1 — F versus (mean) velocity T/{N — 1) for various lengths of 
the chain, 6 = —tt/2. An optimized path x = Xopt{t) (see Fig. U) is used to changed 
the coupling in time. 



5. Experimental implementation 

Let us finally discuss the implementation of the bilinear-biquadratic Heisenberg 
Hamiltonian and the realization of the spatial-temporal control of the couplings. We 
consider atoms with spin F = 1 [i.e., atoms with an 2F+l-dimensional hyperfine degree 
of freedom) at low temperatures confined in a deep (ID) optical lattice. The system 
then is well described via the Bose-Hubbard Hamiltonian [T3l [26] 



Here ai^„ annihilates a particle in a hyperfine state with = a at site i, fii = ^ 



is the number of particles, and Fi = J^acr' Q^lo-^o-o-'^«,o-' is the (total) spin operator on 
site i (T being the usual spin-1 matrices). The tunneling amplitude t is obtained 
from the overlap of the Wannier functions w{x) [26]. Tunneling conserves both, the 
total spin S and the total spin projection m^. The parameters Cq = {go + '^gi)/'i and 
C2 = {g2 — fl'o)/3 depend on the effective ID interaction strengths in the spin 5" channel, 
gs- Those are proportional to the two-body spin-dependent scattering length as- Their 
absolute strength can be controlled via the orthogonal confinement, the relative strength 
moreover can be changed via (magnetic or optical) Feshbach resonances. Note that due 
to the bosonic character of the particles, on-site contact interactions are zero for odd 
total spin. 

The ratios t/co, t/c2 between tunneling and on-site interactions are tunable via the 
lattice parameters (lattice depth and orthogonal confinement). For t <^ |co|, |c2| the 
system is in a Mott-insulating state with atoms being quenched at fixed lattice sites. 
Here we will only consider the case of having a single particle per lattice site [15], and 
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assume tunneling to be sufficiently weak compared to on-site interactions, such that it 
can be treated perturbatively with t/co as small parameter. Then, the following effective 
spin-spin Hamiltonian can be obtained in second order perturbation theory [T3] : 



Co 



1/1 2 - \2' 

JiJi+l ) + 77 + ^ J^Ji+l 



C2 V / 3 \ 1 -|- C2 1 — 2C; 



(20) 



where 62 = C2/C0. This Hamiltonian is equivalent to the one of Eq. ([T]). For ^^Na, the 
bare values of the scattering lengths are very similar: = 46as, 02 = 520^ [27]. Then 
C2 ~ 0.04 and 6 ^ — 0.747r. However, through auxiliary magnetic and electric fields, it 
is possible to modify the system such that the complete anti-ferromagnetic part of the 
phase diagram can be reached [T4] . 

Let us now move to the realization of the adiabatic passage. The necessary spatial- 
temporal variation of the couplings can be realized by two pairs of laser beams of 
wavelengths A and A/2, respectively, and identical polarizations, counter-propagating in 
X direction (additional optical potentials in y- and z-direction are necessary to confine 
the atoms in 3D). The trapping potential in x-direction seen by the atoms then is 
proportional to the square of the electric field. Thus 

Vi^tix, t) oc /half (t) cos2(27rx/(A/2)) + /^[(t) cos2(27ra;/A + </)Mi(t)), (21) 

where Jhaif and Jfuu are the corresponding laser intensities. We have introduced an 
additional phase shift 0fuii for the laser of wavelength A. The effective spin-spin-coupling 
is proportional to t^, see Eq. fl20|) . which in turn has an exponential dependence on the 
height of the potential between adjacent sites. Setting only Jhaif > (and Jfuu = 0), a 
lattice (with a distance of A/4 between adjacent sites) is defined with all couplings being 
equal. We assume this lattice to be loaded with a single particle per site in the ground 
state. Setting </)fuii = n/2, then Jfuu is increased to strongly reduce the coupling between 
each second pair of sites. At the same time, /half has to be adjusted to keep the tunneling 
between the other sites in a regime where the effective Hamiltonian ( l20l) is valid. This 
provides the initial situation for the transport process, c./.. Fig. [6] (a). After preparation 
of the ffist spin, Jfuu is decreased to zero (adjusting properly also Jhaif; see Fig. [6] (b)), 
and 0fuii is set to zero. Then increasing Jfuu again allows to turn off selectively only the 
other subset of spin-spin-couplings [Fig. [6] (c)]. This procedure realizes the change of 
the Hamiltonian necessary for the adiabatic passage transfer process. 

Finally, we can make a rough estimation of the time scales of the process which can be 
achieved with realistic optical lattice parameters. For ^^Na, A/2 = 514 nm, and for tight 
orthogonal confinement, in the perturbative regime (|t/co| ^ 1) couplings a = hJ with 
J ^ 100 . . . 200 s-i could be achieved, leading to Tjv=9 ^ (1 ... 25) ■ lO'^ s for iV = 9 sites 
and an error e < 10~^ in the optimized case. For atoms with spin F = 1, like Sodium, 
deep in the Mott insulating phase, decoherence is caused either by 3-body losses or by 
scattering of photons. Notice, however, that 3-body losses are relevant only for filling 
factors larger than two. If the lasers forming the optical lattice are far detuned from all 
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atomic resonances, the spontaneous scattering rate is highly suppressed and coherence 
times much larger than T/v=9 are easily accessible. That means that the transport scheme 
we propose via adiabatic passage in spin chains should be experimentally realizable. 

6. Conclusions 

In this article we have presented a scheme to transfer a quantum state through a spin 
chain using an adiabatic passage technique. The most remarkable features of this scheme 
are the following. First, the transfer fidelity can be made arbitrarily close to one by 
increasing the transfer time. Second, as long as the change of the couplings is adiabatic, 
the transport fidelity does not depend on the exact path x(t) in parameter space. There 
is thus no need of a precise control of couplings and timing. Third, in contrast to other 
transfer schemes in spin chains, there is no need for a receiver, meaning that once the 
transfer has been accomplished the system remains frozen. 

We have applied this method to a spin-1 chain in the anti-ferromagnetic dimerized phase. 
Using standard optimal control techniques, we have shown that the transfer velocity can 
be substantially increased if an optimized path in parameter space is chosen. We have 
proposed a realistic experimental implementation of adiabatic passage transfer using 
ultracold atoms and optical superlattices. 

In the spin-1 dimerized phase, the transfer velocity is limited by the fact that for a chain 
with an odd number of sites the gap vanishes as ^ oo. Identifying a system with 
similar characteristics but with a gap that persists in the thermodynamic limit would 
allow for finite transfer velocities regardless the size of the chain. 




0.5 1 1.5 2 2.5 



Figure 6. Sketch of the optical superlattice potentials obtained from varying the 
laser intensities /half and If^ii to obtain the required control over the couplings, 
(a) /full > 0, /half > 0, (ptuu = tt/2; (b) /fuu = 0, /half > 0, (/)fuii = tt/2; (c) 
/full > 0, /half > 0, 4>{uu — (see text for details). 
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